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Abstract: In proton-proton collisions at LHC energies, and low mass Higgs bosons 
would be produced with high and predominantly longitudinal boost with respect to the 
beam axis. This note describes a new analysis tool devised to handle this situation in cases 
when such bosons decay to a pair of r-leptons. The tool reconstructs the rest frame of 
the T~^T~ pair by finding the boost that minimises the acollinearity between the visible r 
decay products. In most cases this gives a reasonable approximation to the rest frame of the 
decaying boson. It is shown how the reconstructed rest frame allows for a new method of 
mass estimation. Also a considerable gain in sensitivity to helicity correlations is obtained 
by analysing the r-jets in the reconstructed frame instead of using the laboratory momenta 
and energies, particularly when both r-leptons decay hadronically. 
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1 Introduction 



Experimental precision data strongly suggest that the Higgs boson should be relatively 
light. While experimental data exclude Standard Model Higgs, H, masses below 114.4 
GeV/c^, the lightest neutral supersymmetric Higgs can have a mass as low as 92.8 GeV/c^ [1], 
which is very close to the mass of the In the mass range below 160 GeV/c^ the 
decay mode is of particular importance in searches for neutral H, as demonstrated in Teva- 
tron searches [2]. In this mass region, H production in proton-proton (pp) collisions at 
LHC energies is dominated by gluon-gluon fusion [3], and a search for gg ^ H ^ t^t^ 
production at the LHC is very important, as emphasised in [4]. 

Together these facts strongly motivates a careful study of the t^t~ system as produced 
at the LHC, and the development of analysis tools that enhance sensitivity to measure 
properties of H bosons. Furthermore, since production is an irreducible background 
for this channel it is also important to fully understand the behaviour of Z^^t~^t~ events 
and to look for possible contributions from other processes in a selected sample of r^r^ 
pairs. The sensitivity of such a study can be enhanced by exploiting the fact that the 
parity violating r decays carry helicity information. Thus, r-leptons can be used as spin 
analysers, and provide information that could be helpful in distinguishing the existence of 
a scalar H in a background of decays. 

In [5] some variables that could be used for a spin analysis are proposed. However, 
these are all defined in the rest frame of the decaying boson, a frame that cannot be found 
in any straight forward manner because of the escaping neutrinos produced in the r-decays. 

Methods have been proposed to assign the measured missing transverse energy, ^-p; 
to the neutrinos to obtain unbiased estimates of t~^t~ mass [6, 7]. The resulting estimates 
of the neutrino momenta also gives access to an approximation of the rest frame of the 
decaying boson, but implications of this beyond mass estimation is not discussed. Further, 
these methods only work in a fraction of the events with a suitable topology whereas the 
technique proposed in this paper is applicable for all event topologies. 

This note presents a new simple algorithm that finds an approximation to the rest 
frame of the decaying boson. This will be shown to result in an improved sensitivity to 
spin, and to an alternative variable that could be of use for the estimation of the mass of 
the decaying boson. 

2 Reconstructing the r+r^ rest frame 

Because of their short life time r-leptons are typically studied indirectly through their decay 
products. Due to the one or two neutrinos produced in the decay, the 4-momentum of the 
visible r-decay products, the r-jet, has a broad spectrum, often carrying less that half 
the energy of their mother r. Therefore, reconstructing the rest frame, RF, of a particle 
decaying into a r^r~ system is in principle not possible. 

However, for a sufficiently massive particle decaying into a r+r" pair, the directions 
of the r-jets are expected to be close to that of their mother r-leptons. Thus, one expects 
the r-jet pair to have just only a slight deviation from being back-to-back in the rest frame 
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of the massive particle. It is useful to define the acoUinearity - the angular deviation of 
the r-jets from being back-to-back. This acoUinearity is a well defined positive number in 
any frame of reference. Henceforth, a will denote the acoUinearity between two r-jets, and 
furthermore orf the acoUinearity in the heavy particle rest frame. 

Simulations show that at 7 TeV pp-collisions, the and a light H, are mostly produced 
with a high and predominantly longitudinal boost, (3^. Thus, measured 4-momenta will 
deviate significantly from those in the heavy particle RF. Exploiting that ajip should 
be small and assuming that /?t is small, the method proposed in this paper consists in 
searching for the f3z that minimises a. When restricting /3 to be parallel to the beam axis, 
the a as a function of (3 has a single minimum and hence can be reliably minimised. 

The method can also extended to look for the transverse component of the boost if 
one has a good estimate of the direction of the heavy particle in the transverse plan. The 
transverse direction of the heavy particle, can be estimated by summing up the transverse 
momenta of the r-jets and the reconstructed missing transverse energy, ^-p. Our algo- 
rithm minimises a in two steps; first, by varying the longitudinal boost, /3z, secondly when 
a suitable minimum is found, the transverse boost is varied along the found transverse 
direction. In both steps the a is 1-dimensional function with a single global minimum and 
no secondary minima. 

In this paper, we will denote the two search strategies as the Z- or XYZ-method 
depending on whether we only will try to estimate the longitudinal component of the 
boost of the heavy particle or estimate the transverse component as well. 

Since particle directions are well measure quantities these methods are robust and can 
be applied in every collision event for any pair of measured particles or jets. 

Performance 

The performance of the methods were studied using r+r^ pairs from simulated and 
114 GeV/c'^ gluon-gluon fusion H events produced in 7 TeV pp-collisions. Anticipating 
that backgrounds could be difficult to handle in multi-prong r decays, only decays with 
exactly one charged particle in the final state were considered here. However, the inclusion 
of r decays with more than one charged hadron is not expected to change the performance 
of this method significantly. 10'' Z^ and H events were generated with PythiaS [8] and 
subsequent r-decays were done with TAUOLA [9] to correctly include spin and polarisation 
effects. Detector resolution effects were not included. However, an angular cut requiring 
the two jets to have \rj\ < 2.5 was imposed to match the acceptance of a typical tracker in 
an LHC experiment. 

As shown in figure 1(a) the polar angle distributions for both H and Z^ decays are 
predominantly parallel to the beam axis. Therefore, assuming a boost direction parallel 
to the beam seems to be a reasonable guess for a very large fraction of the events and 
we therefore only expect the XYZ-method to be slightly better in reconstructing the rest 
frame than the Z-method for these events. In figure 1(b) it is seen that the true orf 
is concentrated at small values, with (orf) = 4.5° for Z^ events, while the acoUinearity 
distributions in the reconstructed frames are peaked towards low acollinearities as expected. 
The effect of working in the collinear frame (a = 0) instead of in the true RF should have 
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(a) Polar angle, 6, of the Z° and 114 GeV/c^ 
H bosons. Typically the bosons are produced at 
small angles close to the beam axis. 



(b) AcoUinearity, a, for Z° — >■ r+r^ in the gener- 
ated and reconstructed RFs and in the laboratory 
frame of the Z'K 



Figure 1. Angular distributions for simulated heavy boson decays. 



a small effect if the direction of the boost is correctly estimated. In the collinear frame 
closest to the true RF, the momenta of the two r-jets will deviate with a factor cos(aRF) 
from the truth, meaning a typical deviation of 0.3%. 
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Figure 2. Average acoUinearity as a function of attempted longitudinal boost for various ranges 
of true boost. 



Since the a is always positive, it is easy to minimise with respect to the applied /3z_t 
of the two tau jets. For a selection of ranges in the generated /3 of a produced Z'^-boson, 
figure 2 shows how, on average, a varies as a function of an applied longitudinal boost to 
the two jets. The minimal acoUinearity is always at a longitudinal boost that is close to 
the generated one. The actual value of a is reflecting to which extent there might be a 
transverse boost present, ignored in figure 2. 

To estimate the goodness of the reconstruction of the RF, the vectorial difference 
|/3gen — Pvecol is shown in figure 3(a) for both reconstruction methods. The reconstructed 
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(a) Vectorial difference between the generated and 
reconstructed RF. 
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(b) Difference in magnitude between the gener- 
ated and reconstructed /3. 



Figure 3. Comparisons of both reconstructed and generated rest frames in Z'^ t^t events. 

RFs are seen to be closer to the generated RF for the XYZ-method as expected. Figure 3(b) 
shows the difference in magnitude between the generated and reconstructed boost with a 
FWHM of the distribution of 0.010 for the XYZ-method. 

In conclusion both methods reconstructs the boost very reasonably, with a distri- 
bution around the true value with a mean deviation around 0.076 for the XYZ-method 
and 0.149 for the Z-method, when averaging over all events. Due to the exclusion of the 
transverse component, a small offset away from zero in the most probable value is seen 
in the Z-method. Although figure 1(a) shows that the bosons are mostly boosted longi- 
tudinally, improvements are found when also estimating the transverse components of the 
boost. However, these improvements will be very dependent on the ability to estimate the 
transverse direction of the bosons correctly in the detector. Therefore the two methods 
act complementary; for events with a small transverse boost and low precision on the 
the Z-method should be applied whereas for events with good precision in ^-nd a high 
transverse boost improvements will be gained by using the XYZ-method. 

3 Estimating the Boson Mass 

Due to the neutrinos in the r decays, the visible invariant mass, Mvis, distribution peaks 
far below the real mass of the resonance produced. A method correcting for this is the 
so called collinear approximation [6]. This method assumes the neutrinos are collinear 
with the decaying r-lepton, and includes .^x by projecting it to the r-jets. However, this 
projection is only possible and reliable in a fraction of events with large transverse boost 
and well aligned [10, 11]. 

In [7] the performance of a likelihood based technique is presented for simulations of 
Higgs and events in antiproton-proton collisions at 1.96 TeV, and also for data collected 
by the CDF experiment. For fully hadronic decaying tati-leptons, unsmeared simulations 
show a mass reconstruction resolution of 8%, and more realistic simulation shows a mass 
resolution of 14%, and finally the CDF — t- t~^t~ events show a resolution of 16% with 



-5- 



this technique. The smeared results are obtained for events within \r]\ < 1, smearing of 
hadron momenta of 10% and of of 5 GeV in each of the transverse components. 

At the LHC the events are much more boosted, and the experiments have larger r] 
coverage. In [12], CMS expects a likelihood mass resolution of 21% for a 130 GeV Higgs 
(for leptonic and semi-leptonic r-decays), in an analysis covering \ri\ < 2.3 for hadronic r 
decays. This is approaching the width of the the Mvis distribution which is reported to be 
24%. A possible explanation for the poorer resolution in CMS could be that rneasures 
a smaller fraction of the full neutrino energies for r-jets with large tj than for transverse 
ones. Hence, estimates of the full neutrino energies from ^-p arises from correspondingly 
larger scale factors. 

Another mass variable for r pairs is the effective mass, Mcfj = ^ {Pr^ ^ Pt'^ 
where Pr are the visible r-jet 4-momentum, and is the missing transverse momentum 
vector. In addition to the M^is and Mefj recent result from the ATLAS collaboration [13] 
has also used a likelihood based method, whereas the recent result from the CMS collabo- 
ration [14] simply uses Mvis in the analysis. 

In the following an alternative variable for mass estimation is described and for sim- 
plicity comparisons have been made to distributions of Mvis which is a straightforward 
robust alternative in use by all experiments. 

Finding the kinematic edge 

The mass of a particle decaying to r+r^ could be inferred by determining the kinematic 
endpoint of the distribution of twice the leading r-jet energy in the boson RF. This quantity 
will henceforth be named as 2E^-p- opposed to the collinear approximation this quantity 
can be calculated for all events. In figure 4 distributions of 2£'j!^p are compared to Mvis 
distributions. 

This way of estimating the mass becomes particularly useful when both r-leptons decay 
into hadrons, where only one neutrino is present in each r decay. Figure 4 suggests that 
a first guess to the shape of 2£'p^p could be a triangle convoluted with some resolution 
function, although it is clear that the inclusion of effects like the / H width and helicity 
effects modifies such an assumption. Still the kinematic edge can be found be determining 
the point of steepest slopes of the distribution. 

Figure 4 is the result for all hadronic decaying events without any selection. As 
mentioned earlier, the decision on which boost reconstruction method to use should depend 
on the magnitude of the transverse boost. Therefore, the events have been divided into 
two classes; one with a large difference in azimuthal angle between the two jets, Ac/^)®**^, 
compatible with no transverse boost, where the Z-method were used, and one with the 
rest of the event where the XYZ-method were applied. Here a threshold of lS.(j)^'^^^ = 170° 
were chosen to separate the two classes. The criterion for the choice of boosting method 
would have to be subject to study in a true experimental environment. Furthermore, a cut 
of 10 GeV on the for both r-jets were applied to both classes. In figure 5, it is shown 
that 24% of the Z'^-boson and 16% of the Higgs events can be reasonably subjected to the 
Z-method. Hence, without any use of approximation of the centre of mass system 

is still available of these events. 



-6- 



20 40 60 80 100 120 140 160 180 200 

[GeV] 



20 40 60 80 100 120 140 160 180 200 

[GeV] 



(a) Distributions for events. (b) Distributions for H events. 

Figure 4. Distributions for ^E^^ calculated using the XYZ-method and invariant mass of the 
visible r-jcts, Mvis, for t^t~ pairs both decaying to single charged hadrons. The dotted line marks 
the position of the steepest slope for the 2i?^p distribution. 
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(a) Class 1: Distributions for events with (b) Class 2: Distributions for events with 

Figure 5. Distributions for 'i.E^^ calculated for the two event classes using the Z and XYZ-method 
for r+r" pairs both decaying to single charged hadrons. 



Since the distributions of 2E^y ^'^^ clearly asymmetric and non-gaussian, the spread 
of the distributions calculated from the most likely value are therefore also asymmetric. 
For Higgs events with /S.(j)^'^^^ above 170°, the spread from the peak position is found to be 
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To evaluate the precision of retrieving the boson mass by finding the kinematic edge, 
-^boost) pseudo-experiments were performed. From the 2£'p^p distribution shown in fig- 
ure 4(a) several pseudo-distributions were generated with varying number of events and 
for each distribution the kinematic edge was found as the point of the steepest decent. 
In table 1, the found values and spread of Mboost ^-re listed as a function of the number 
of events in the pseudo-experiment for the two classes of events, where class 1 refers to 
the events in which the Z-method has been used, and class 2, refers to the events where 
the XYZ-method were used. For the simple kinematic edge-finder implemented for this 
analysis the Mboost converges to 90.8 it 1.4 GeV for 4000 events (for class 1 events). These 
values indicates that the method could be of interest to investigate further in experimental 
data. 

For a triangle function convoluted with a gaussian, numerical calculations have shown 
that taking the kinematic edge to be the point of the steepest slope overestimates the 
edge value. For a gaussian width of 7-10 GeV, the systematic shift is found to be 0.5- 
1.0 GeV. For an experimental analysis including backgrounds a more sophisticated edge 
finding method, e.g. template fitting, should be deployed. 

Table 1. Mass estimation A/boost as a function of number of events ricvt as obtained from pseudo- 
experiments for the two classes of events. Class 1 consists of events with A0j°*^ > 170° where the 
Z-method have been applied and class 2 consists of events that has been subject to the XYZ-method 





Obtained values of Mboost 


Number of events 


Class 1 


Class 2 


250 


88.3 ± 5.3 


88.5 ± 4.7 


500 


89.3 ± 3.1 


89.3 ± 2.5 


1000 


90.3 ± 2.3 


89.7 ± 1.7 


2000 


90.6 ± 1.8 


89.9 ± 1.3 


4000 


90.8 ± 1.4 


90.0 ±1.1 
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Figure 6. Correlations between Mvis and ^E^p using the XYZ-method in t^t" events with 

all decay channels of the r-leptons included. 
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Intuitively, constructing a mass variable that accumulates towards the real mass should 
help distinguishing different masses. One way to distinguish distributions is to compare the 
largest difference of the cumulative distributions, which is input to the Kolmogorov- 
Smirnov test. It is found that while between the calculated ^E^^ distributions (using 
the XYZ-method) for and H events is 0.35, it is only 0.24 when comparing the two 
Mvis distributions. It is not clear if the tails towards high values of 2E^p seen to be more 
important in figure 4 than figure 4(a) is helpful in a search scenario, because clearly the 
contribution is non-negligible for values well beyond the kinematic edge. At present, it 
cannot be concluded that this variable will be better than any of the other proposition in 
a real experiment. However it seems to be worthwhile to consider the use this variable in 
an experimental search for the Higgs boson. 

In figure 6, the correlation between Mvis and 2E^-p is shown for Z^ events. Thus, since 
the correlation is significantly lower than 1, the 2i?p^p variable is adding information, and 
it should be studied for use in a search scenario. 

4 Helicity correlations 

Several powerful variables for studying r polarisation in Z^ decays were developed for the 
LEP experiments [15]. These variables are defined in the rest frame of the decaying boson. 
While this is close to identical to the laboratory frame at LEP, this is not true at the LHC, 
and this paper proposes to study these variables in a reconstructed RF instead. 

For —7- vr^ Ur events, the distribution of the decay angle in the r-lepton rest frame, 
that is the angle between the direction of the r-lepton and its visible decay product, is 
determined by the r helicity. In the rest frame of the heavy boson producing a t~^t~ pair 
this translates to a well determined distribution in the fraction of r energy carried by vr. 
The quantity used in the following figures is the energy of visible particles in the decay, 
i?™ , normalised to the kinematic endpoint (half the mass of the decaying boson) . Weaker 
energy-energy correlations also shows for the — t- 1/^ vi and — t- /i^ Vt nh^ modes. 




Figure 7. 2E^^^/AIz for r-leptons with positive and negative helicities decaying into n^v-r shown 
in the generated and reconstructed RF using the XYZ-method as well as the laboratory frame. 
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The E^^'^ distributions for single vr decay modes shown in figure 7 for the generated 
and reconstructed heavy boson RFs show that the true correlations are partly recovered 
in the reconstructed RF using the XYZ-method proposed in section 2. Furthermore, it is 
shown that these correlations are much weakened in the laboratory frame. 

4.1 Analysis of the sensitivity to spin 

Conservation of angular momentum leads to distinctly different helicity configurations when 
comparing the final states of and H decays in their rest frames. While the spin compo- 
nents along the direction of flight of the r"*" must add up to ±1 for the Z^, the sum must 
be for the H decays. Thus the spin effects can be studied by looking at the correlations 
in i?™ of the two r-leptons as shown in figure 8. 




2E"/M^ 2E"/Mj 



(a) Spin (H sample). (b) Spin 1 (Z" sample). 

Figure 8. Energy correlations in the reconstructed RF using the XYZ-method of a r+r" pair with 
both T-leptons decaying to Vt . To avoid effects coming from the mass differences all energies 
are scaled with the mass of the decaying boson. 



To quantify the observed difference between spin and spin 1 particles, the Z and H 
samples were split into training and test samples. Probability density functions for both 
spins, Po{E^- , E^+) and Pi{E^~ , E^+) , were constructed as the fraction of events in a bin 
of {E^^,E^+) in the corresponding training sample. From this a likelihood, £, is created 
for N events from both training samples by 

N 

log £ = ^ (log(P{) - log(P^)) /N (4.1) 

i=0 

The test samples of Z^ and H events are divided in sub-samples of events and the 
likelihood is calculated for each sub-sample. The distribution of logC when N = 1000 and 
when both r-lepton decays to Vt is shown in figure 9. 

Inspired by [15], a sensitivity^ S, is defined such that the number of standard deviations 
between the means of the two distributions, n^, is given by 

n„ = sVn (4.2) 
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Figure 9. Log likelihood plot for energy-energy correlations in the reconstructed RF using the 
XYZ- method when both r-leptons decays into Vt- 



Table 2 shows values obtained for S in the reconstructed RFs using the Z- and XYZ-method 
as well as in the laboratory frame for the different r decay modes. For all decay channels 
sensitivity is gained when using the reconstructed RF over the laboratory frame. When 
just applying the Z-method, differences in the transverse momentum of the bosons are 
kept in the distributions, while the XYZ-boost approximates the RF more correctly. The 
Z-method thus analyses a combination of the two effects, the boson and spin, while the 
sensitivities determined using the XYZ-method are largely only due to differences in spin. 



Table 2. Sensitivity to spin of the boson for different decay modes of the two r-leptons shown in 
the true and reconstructed RFs as well as in the laboratory frame. 



T decay 


mode 


Sensitivity 


n 




Reco.RF(Z) 


Reco.RF(XYZ) 


Lab. 


True RF 


Vr Vl 


Vr Vl 


0.13 


0.09 


0.06 


0.08 


£^ Vr ve 


Vr 


0.17 


0.13 


0.06 


0.15 


£^ Vr Ve 


Vr nh^ 


0.15 


0.07 


0.05 


0.13 


Vr 


Vr 


0.42 


0.47 


0.23 


0.51 


Vr 


Vr nh^ 


0.21 


0.15 


0.10 


0.18 


Vr nhP 


Vr nh^ 


0.19 


0.13 


0.09 


0.18 



Finally, further sensitivity would be gained by using ratios between neutral and charged 
energy for Vr nh^ final states. An exhaustive analysis of spin should probably make 
use of these ratios. These ratios remain almost unchanged in any boosted frame. The 
purpose of this note is to discuss the advantages of finding an approximate rest frame, so 
the sensitivities when using these variable are not reported here. 
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5 Conclusions and discussion 

By minimising the acollinearity, two methods to reconstruct the rest frame of boosted heavy 
particles decaying to two r-leptons have been proposed, and it has been demonstrated that 
they both find the rest frame reasonably well. Whereas the XYZ-method requires an 
estimate of the transverse direction of the heavy particle, the Z-method can be used for 
analyses where one does not have a good estimate of this direction. 

A new technique for mass estimation has been proposed using the end-point of the 
leading r-jet energy in the reconstructed rest frame. This technique does not introduce 
assumptions on transfer functions or production mechanisms of the heavy particle as more 
elaborate reconstruction techniques, and it can be applied to all events within the geomet- 
rical acceptance of the experiments. The method is complementary to existing methods in 
the sense that it can be applied without using information on the measured -^Xi £^^d since 
it is not directly correlated with the visible mass it adds information, and thus the two can 
be combined to gain further sensitivity. 

Furthermore, a scheme for an event-wise spin analysis of the t~^t~ system has been 
outlined, and it has been shown that enhanced sensitivity to spin is gained when trans- 
forming the 4-momenta as proposed in this paper. Both boost methods presented enhances 
sensitivity to spin when compared to an analysis in the laboratory frame. Which method 
to use depends on whether or not differences in production mechanism of the boson should 
be taken into account. 

In summary, a tool has been developed to approximate the rest frame of t~^t~ systems, 
and two applications thereof have been illustrated. These applications could prove to be 
helpful in the analyses of decays or in searches for new particles like H, Z' or SUSY 
particles. Additionally, it is possible that a transformation to the rest frame could be 
beneficial to other variables used in the study of r^r^ systems. 
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